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Abstract 

There are several numerical methods for computing approximate zeros of a given univariate polynomial. In 
this paper, we develop a simple and novel method for determining sharp upper bounds on errors in approximate 
zeros of a given polynomial using Rouche's theorem from complex analysis. We compute the error bounds 
using non-linear optimization. Our bounds are scalable in the sense that we compute sharper error bounds for 
better approximations of zeros. We use high precision computations using the LEDA/real floating-point filter 
for computing our bounds robustly. 
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1 Introduction 

The problem of solving the polynomial equation 

g(z) = ao + a\z + d2Z 2 + ... + a n z n = (1) 

is a fundamental problem in algebra. Here, the coefficients m, < i < n are complex numbers. It is now folklore 
that no closed form formula with arithmetic operations and radicals exists for computing zeros of polynomials 
with degree n > 4. In spite of the absence of solution formulae, the fundamental theorem of algebra states 
that equation ^ always has n zeros, real or complex. The subject of computing or approximating these zeros has 
been called algebraic aspects of the fundamental theorem of algebra ^21- The computational challenge is in 
determining zeros precisely for high degree univariate polynomials. In computer algebra applications, one usually 
needs to solve equation for large values of n, typically well above 100 or of order of several thousands Jlj. High 
precision of hundreds of bits is frequently required for the representation of the coefficients a , ai, a n and the 
zeros. Development of robust and efficient solutions for equation with large n requires much research. 

1.1 Computing approximate zeros of polynomials 

Since there is no closed form formula for finding zeros of a polynomial of degree greater than four, we need to 
find zeros by numerical methods. The first iterative algorithms with guaranteed convergence to all the n zeros 
of g(z) (for any input polynomial g{z) of degree n ), are due to Brouwer and Weyl, both published in 1924 (see 
[l4*|b Pan's new 0(n 2 lognlog&n) algorithm to approximate all the zeros of a n-degree polynomial within an error 
bound of 2~ b maxj\zj\, where Zj are the zeros of the polynomial is the best known result so far. There are 
several algorithms for finding zeros of polynomials based on Newton's method Laguerre's method ^Hl> Jenkins 
and Traub's method and Muller's method [14 . 

1.2 Estimating upper bounds on errors in approximate zeros of a given polynomial 

Since all the above methods essentially compute only approximate zeros, we need to know how close these approx- 
imations are to the actual zeros. Consider a practical scenario where we have pretty good estimates of all the zeros 
of a known univariate polynomial as in equation ^ Although the zeros may be quite accurate, we still need to 
know how good they are. In other words, we need to know how much error is present in each approximate zero. 
Consider the problem of developing a method for determining upper bounds on the errors in each approximate zero, 
irrespective of the method and computational precision by which such approximate zeros have been determined 
apriori. Smith jTH] computes upper bounds on errors in approximate zeros using a method based on Gerschgorin's 
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in given approximate zeros of a given univariate polynomial. Our method computes an upper bound on the dis- 
tance of an approximate zero from the exact zero given all the coefficients of the univariate polynomial and all the 
approximate zeros. Note that we are given only the univariate polynomial and its approximate zeros as inputs; we 
may not know any actual zero. Towards this end, we use Rouche's theorem from complex analysis [2] as the basis 
for our approximation scheme. Based on this scheme, we develop a non-linear optimization step to determine a 
sharp upper bound on the modulus of the distance of an approximate zero from the (possibly unknown) exact zero. 
This method is reported in |151 flfi] . The implementation of our method for estimating these error bounds uses the 
LED A/real floating-point filter 0^3] so that computations can be done at requisite high precision beyond double 
precision, yielding correct results in the evaluation of inequalities involving arithmetic expressions. The number of 
significant digits in the approximate zeros and the coefficients of the given polynomial are not restricted therefore. 

1.3 Preliminaries 

We need the definition of an analytic function [2. . A complex function f(z) is said to be analytic in a region if it is 
defined and differentiable at every point z in the region. Now we state Rouche's theorem. 

Theorem 1.1 j2j0] Suppose the functions f(z) and g(z) are analytic inside and on a simple closed curve C. If / 
and g have no zeros on C and \ f(z) — g(z)\ < \f(z)\ for all z on C, then the functions f(z) and g(z) have the same 
number of zeros inside C. 

We know that all polynomials are analytic in the complex plane. We know from the fundamental theorem of algebra 
that any polynomial with complex coefficients of degree n has n (possibly) complex zeros. Given a polynomial of 
order n and its n approximate zeros, we wish to bound the error in each of the n aproximate zeros. We present 
our method in the next section. 



2 A novel posteriori error bound for an approximate zero 

Let the given polynomial be 

g(z) — a + a\Z + a 2 z 2 + ... + a n z n (2) 

where the coefficeints a^, < i < n are complex in general. Let the approximate zeros of g{z) (say, as computed 
by some numerical method) be ai, oli, ■ ■ ■ , a n . Let 

n 

/0) ^a n \{{z- ai ) (3) 

i=l 

be the approximation of g(z) as constructed from the given approximate zeros. Now consider the error polynomial 
h{z) = f{z)-g{z). Let 

h(z) = b + b 1 z + + 6„_iz™ 1 (4) 

Note that the error polynomial h{z) has degree at most n — 1, less than that of g{z). Now we apply Rouche's 
theorem on the polynomials f(z) and g(z). Consider ctj, an approximate zero. In this section and in sectional we 
concentrate on ctj] the arguments in these sections apply to each approximate zero. We draw a circle C with a 
small radius r around aj. If the error polynomial h(z) = f(z) — g(z) and the computed polynomial f(z) satisfy 
Rouche's condition \ f(z) — g{z)\ < \f(z)\ for all z on the on this circle, then f(z) and g(z) have precisely the same 
number of zeros inside this circle. If there is no other approximate zero ai, i ^ j within that circle, then f(z) and 
g(z) have precisely one zero in the circle. So, the unique actual zero of g{z) is approximated by ctj and is within 
that circle of small radius r with centre at ctj . Here r is an upper bound on the error in the approximate zero aj . 

For enforcing the condition in Rouche's theorem we observe that, for any z on the circle C of radius r, \z\ — 
\z — aj + aj\ < \z — aj\ + \atj\ = T + \aj\ = R, say. From equation^) \h(z)\ < |6 n _i||z™ _1 | • ■ ■ + |6q|- Hence on C, 

\h(z)\<\b n - 1 \R n - 1 + --- + \b \ (5) 

Also, for all z on the circle C of radius r, \z — cti\ = \(z — aj) — (ai — <Xj)| > \\z — aj\ — \ai — a^- 1 1 = \r — \ai — aj\\, 
Also, \f(z)\ = \a n \\z — a\\\z — a 2 \ ■ ■ ■ \z — a n \, from equation^ Hence on C, 

n 

|/(*)|>|on|r 11 \r-\t*i-aj\\ (6) 

Finally, from inequalitiesElandEl we observe that |o n |r| n"=i i^j \ r ~\ a i~ a j\\ > + + \bo\ implies 

Rouche's condition \f(z) — g(z)\ < \f{z)\ for all z on the entire circle C. We can write this condition as 

l(r) 

r >--+-(= q(r), say) (7) 
m(r) 
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degree n — 1 and that g(r) is a meromorphic function. 

Now we summarize the main theorem of this paper as follows. 

Theorem 2.1 Let g(z) be a polynomial of degree n with n distinct zeros. Let a\, ot2, ■ ■ ■ ,a n be approximations 
to the n zeros of g(z) and let f(z) = a n Y[i=i( z ~ a i) the polynomial with zeros a\, ot2, • • • , a n . Then, for any 
1 < i < n, the error in the jth approximate zero otj is bounded by any real number r satisfying the inequality 
r > q{r), provided \a>j — cti\ > r, for all 1 < i < n where i ^ j. 

Proof. As shown above, Rouche's condition is satisfied on the circle C centred at ctj with radius r provided 
r > q(r). If no other approximate zero ai, i ^ j is in the interior of this circle C then f{z) has a single zero ay 
in the interior of C. By Theorem ll.il g{z) too has a single zero (say, Zj) inside the circle C, yielding the upper 
bound r on the error \aj — Zj\, 



3 A method for computing error bounds 

We state the inequality as 



r * m(r) 

We observe that the inequality (and therefore the inequality 7) is indeed satisfied for sufficiently large values 
of r; this is due to the fact that r * m(r) is a degree n polynomial and f(r) is degree n — 1 polynomial. We state 
this a follows. 

Lemma 3.1 There exists a value or r satisfying inequality |H1 

However, we are interested in finding very small values of r satisfying r > q(r). In particular, note that in our 
Theorem 12 .11 such an r is an upper bound on the error in an approximate zero ctj provided \atj — <Xi\ > r, for all 
i ^ j. We proceed to develop a method to find such a value r as follows. 



3.1 Computing the error bounds: Algorithm I 

We argue that inequality is satisfied for values of r close to q{0), provided f(z) is a close approximation of g(z) 
and the given approximate zeros are well separated. If f{z) is close to g(z), then 1(0) is small and therefore q(0) 
could be small. By Taylor's expansion, q(r) = q(0) + rq'(O) + ^r(?"(0) + • • •. Since higher derivatives of q(r) will be 
small for well separated approximate zeros, we can neglect higher order terms and write the fixpoint of the function 
q as r = q(r) = q(0) + rq'(0) => r = j~7^ ■ Note that q'(0) is also small since q(r) is a function with a numerator 
l(r) with very small coefficients, and a denominator m(r) with large coefficients. Therefore, the fixpoint of q lies in 
the vicinity of q(0), provided the approximate zeros are well separated and f(z) is a good approximation of g(z). 

Example 1 

g(z) = 100000s 4 + 305000s 3 + 410100s 2 + 310205s + 105105 



Comparison with Smith's bounds 



Actual 


Zeros by 
ZERPOL 


Smith's bounds 
on ZERPOL 


Value 
of g(0) 


Our bounds on 
ZERPOL zeros 
using Algorithm I 

c = 0.00001 


-1.05 


-1.0500001610 


6.44S - 07 


1.116840490 
680S63B - 06 


1.11687399 
61584E - 06 


-1 


-0.9999998510 


5.97E - 07 


1.08266830 
781336.E - 06 


1.08270078 
81847S - 06 


-0.5 + iVO.751 


-0.5 + iQ. 86660 
2562368 


2.55E - 08 


3.04704508 
007204E - 08 


3.04707641 
434807B - 08 


-0.5 - i\/0.751 


-0.5 - 10. 86660 
2562368 


2.55E - 08 


3.04704594 
388862B - 08 


3.04707641 

434806.E - 08 



Bounds using Algorithm I for e= 0.0001 



Range zeros at 
precision 7 


Value of (?(0) 


Our bounds on Range 
zeros at precision 7 


Number of 


-1.05 


5.42072490014779E - 06 


S^lSO^T'.H^n^-E - 06 


2 


-1.000000 


5.43336060444733E - 06 


5.43444711286507.E - 06 


2 


-0.5 + iO. 8666026 


1.5286312353107E - 07 


1.528781.0!)S1.:M.2:5.E - 07 


1 


-0.5 - iO. 8666026 


1.52863125747657E - 07 


1.52878412060232E - 07 


1 



Bounds using Algorithm II 



Range zeros at 
precision 7 


Starting 
Value 


NR bounds on Range 


Number of iterations 
NR + Algorithm I 


-1.05 


0.01 


5.42194046573675.E - 06 


4 + 1 


-1.000000 


0.001 


5.434582770616771? - 06 


4 + 1 


-0.5 + iO. 8666026 


100 


1.52878500883089-E - 07 


4 + 1 


-0.5 - iO. 8666026 


IE - 05 


1.52878r,00S>iriNi)E - 07 


3 + 1 
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I Zeros at precision lb | 



Bounds using Algorithm I for c= 0.00000001 



Range zeros at 
precision 16 


Value of q(0) 


Our bounds on Range 
zeros at precision 16 


Number of 
iterations 


-1.05 


4.64160094696633E - 15 


4.64160099338234E - 15 


1 


- 1 .000000000000000 


4.51503828651945E - 15 


4.51503833166985E - 15 


1 


-0.5 + iO. 86660256 173 1732 


1.2706ri22!> Lr,72H»E - 16 


1.27065230727944E - 16 


1 


-0.5 - iO. 866602561731732 


1.2706r>22!t ir ! 72!.)2£' - 16 


1.27065230727944E - 16 


1 



Bounds using Algorithm II 



Range zeros at 
precision 16 


Starting 
Value 


NR bounds on Range 


Number of iterations 
NR + Algorithm I 


-1.05 


IE - 10 


4.64160099338286E - 15 


2 + 1 


-1.000000000000000 


0.0001 


4.51503833167034E - 15 


3 + 1 


-0.5 + iO. 86660256 173 1732 


IE - 22 


1.27065230727944B - 16 


2 + 1 


-0.5 - iO. 866602561731732 


IE - 14 


1.27065230727944B - 16 


2 + 1 



Example 2 

g(z) = 1000z 10 - 2500z 9 - 4608002 8 - 9133400z 7 - 50761800z 6 - 88653100* 5 - 
53510400* 4 - 37313000z 3 - 197170000* 2 - 364800000z - 198000000 



Comparision with Smith's bounds 



Actual zeros 


Zeros by 
ZERPOL 


Smith's bounds 
on ZERPOL 


Value of 


Our bounds on 
ZERPOL zeros 
using Algorithm I 
e = 0.00001 


30 


30 


3.34B - 13 


6.16818594 
879537E - 11 


6.1707644 
473737E - 11 


-10 + ilO 


-10 + ilO 


8.3E - 13 


3.99156437 
089776E - 10 


3.99160428 
654147E - 10 


-10 - HQ 


-10 - £10 


8.3E - 13 


3.99156437 
089776S - 10 


3.9916042 
8654147E - 10 


-5 


-5 


2.91E - 10 


3.18914558 
60219E - 09 


3.18914637 
774777E - 09 


l + i 


1 + i 


6.81E - 14 


3.72504031 
303826E - 10 


3.7250775 
6344139E - 10 


1 - i 


1 - % 


6.81E - 14 


3.72504031 
303826S - 10 


3.7250775 
6344139E - 10 


-1 + iVl-2 


-1 + il-0954 
451114 


3.62E - 08 


2.9138084 
647204E - 09 


2.91383760 
280512E - 09 


-1 - iVl.2 


-1 - £1.0954 
451114 


3.62E - 08 


2.9138084 
647204E - 09 


2.9138376 
0280512E - 09 


-1.5 


-1.5 


1.68E - 13 


1.17048091 

956553S - 08 


1.1704926 
2437473E - 08 


— 1 


— 1 


8.25E - 14 


9390312S - 09 


8273006E - 09 
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I zeros at precision r | 
Bounds fr«:>m Algorithm I for c = 0.0001 



Range zeros 
at precision 7 


Value of 

q( D ) 


Our bounds on Range 
zeros at precision 7 


Number of 
iterations 


30 


7 . 2 4 3 1 ; S 7 9 7 7 - j - j 7 E — 09 


7 . 024438!.t8!.t ['.)'>{):> E — 09 


1 


— 1 + i 1 


4 . 543764 492072 -19 E — 08 


4 . 5438099297 1 74 1 E — 08 


1 


— 10 - j 10 


4. 54370 1 1020 72 1 9 E — 08 


4 . 5438099297 1 74 1 E — 08 


1 


— 5 


3 . 622 1 00330 12892 E — 07 


3 . 622 142 5 5 7 192 2'JE - 07 


1 


1 + i 


4 . 2 4 3 ( i 9 1 - j 1 - • 1 9 8 E — 08 


4 . 2404 1 1 5 .)80565 2 E — 08 


1 


l — i 


4 . 2 4 3 9 1 5 4 3 4 9 8 E — 08 


4.240111 5 5 .8 5 5 2 .E — 08 


1 


— 1 + il. 095445 


3 . 3 169 1 0559! ) 7502 E — 07 


3 . 3 1 ( .1 '■) 1 3 7 2 8 1 8 1 2 1 E 1 1 7 


1 


— 1 — i 1 . 9 5 4 4 5 


3 . 3 169 1 055907562 E —■ 07 


3.31694372818121E — 07 


i 










-1 


9.26483214486361E - 07 


9.26492 17931*505 - 07 


1 



Bounds from Algorithm II 



Range zeros at 
precision 7 


Starting 
value 


NR bounds on Range 
zeros at precsion 7 


Number of 
NR +Algorithm I 


30 


IE — 5 


7.02507120754467E - 09 


3 + 1 


-10 + ilO 


IE — 1 


4.54421907500478.E - 08 


4 + 1 


-10 -- ilO 


IB - 10 


4.54421907500478B - 08 


3 + 1 


— 5 


IE - 1 


3.62247181860654.E - 07 


4 + 1 


1 + i 


IE — 5 


4.24079393150912.E - 08 


3 + 1 


1 - i 


IE — 5 


4.24079393150912.E - OS 


3 + 1 


-1 + il. 095445 


IE — 10 


3.317217877106.E - 07 


3 + 1 


-1 - il. 095445 


IE — 2 


3.317217877106.E - 07 


3 + 1 


-1.5 


IB - 2 


1.33255239212745B - 06 


4 + 1 


-1 


IB - 2 


9.26581422929994B - 07 


4 + 1 



| Zeros at precision 16~| 
Bounds from Algorithm I for c = 0.00000001 



Range zeros at 
precision 16 


Value of 

9(0) 


Our bounds on Range 
zeros at precision 16 


Number of 


30 


1.8329<}:;7~.2^~.-.0E .1 7 


1.83299377121549E - 17 


1 


-10 + ilO 


1.18568546630189E - 16 


1.18568547815876E: - 16 


1 


-10 - ilO 


1.18568546630189S - 16 


1.18568547815876E: - 16 


1 




9.41580700240742E - 16 


9.451807u;Hi!.t8.-,r,3-E - 16 


1 


1 + i 


1.10651507753055B - 16 


1.10651508859572B - 16 


1 


1 - i 


1.10651507753055E - 16 


1.1065ir,ii88r J !)r J 72£' - 16 


1 


-1 + il. 095445 
115010332 


8.65540430075807E - 16 


8.65540438731223.E - 16 


1 


-1 - il. 095445 
115010332 


8.65540430075807S - 16 


8.65540438731223fi - 16 


1 


-1.5 


3.47688796630911E - 15 


3.47688800107798E - 15 


1 


-1 


2.41763729896059S - 15 


2.41763732313696.E - 15 


1 



Bounds from Algorithm II 



Range zeros at 
precision 16 


Starting 
value 


NR bounds on Range 
zeros at precision 16 


Number of 
NR. +Algorithm I 


30 


IE - 14 


1.83299377120949B - 17 


2 + 1 


-10 + ilO 


1000 


1.1850854781587E - 16 


3 + 1 


-10 - 1 10 


IE - 04 


1.18568517fsl5fs7.E - 16 


3 + 1 




IE - 04 


9.45180769698555S - 16 


3 + 1 


1 + i 


IE - 14 


1.10651508859572E - 16 


2 + 1 


1 - i 


IE - 14 


1.10651508859572S - 16 


2 + 1 


-1 + il. 095445115010332 


IE - 14 


8.055 10 13fs731227E - 16 


2 + 1 


-1 - il. 095445115010332 


IE - 04 


8.65540438731227E - 16 


3 + 1 


-1.5 


IE - 03 


3.47688800107808E - 15 


3 + 1 


-1 


IE - 03 


2.41763732313701E - 15 


3 + 1 



So, we start testing r with an initial value q(0). We increase r until r > q(r). Note that q(0) > 0, and therefore, 
the inequality r > q(r) is not satisfied at r = 0. We use q(0) as a starting point for searching a value of r satisfying 
r > q(r). In the search, we multiply the value of r by a factor of 1 + e repeatedly (for a small e), until the inequality 
r > q(r) is satisfied. Further, if r < \a>j — Oi\ for all 1 < i ^ j < n, we can assert using Theorem 1241 that r is an 
upper bound on the error in the approximate zero otj . 

We claim that the number of times q(r) is evaluated until we get a solution 5 satisfying S > q(S) is 0(log 2 S). 
This follows from the fact that we start with an initial approximation q(0) and multiply r by a factor (1 + e) until 
the inequality is satisfied, as guaranteed by Lemma Ei.ll We state this as a theorem. 

Theorem 3.2 The value of r (say S), satisfying the inequality can be determined in log 1+e ((5/q(0)) evaluations 
of the function q(r). 

Each evaluation of q{r) involves evaluation of two (n — l)-degree polynomials, l(r) and m(r). Note in the tables 
(Tables 1 through 6) that as long as we select a good starting value q(0) for r, we practically need to execute just 
one step to get to a value of r satisfying r > q(r). Indeed, the error bounds are just a little bit higher than q(0) 
in all cases, very much as argued above. So, we observe that for good approximations of well separated zeros (as 
considered in all our examples), Algorithm I converges very fast. Algorithm I requires larger number of steps if (i) 
the initial value of r is far from a feasible solution for inequality 7 and, (ii) for smaller values of e. 

3.2 Using requisite high precision for computing the error bounds 

Another factor in the time complexity of Algorithm I is the precision at which computations would require to be 
performed to correctly determine satisfaction of inequality 7. We use the LEDA/real floating-point filter (3l ITo|. 
for computations at requisite higher precisions. LEDA first tries to check the inequality using double precision. If 
a decision can not be made at double precision, LEDA uses higher precision to check the inequality. Precision is 
increased until a decision is possible. High precision computing is necessary since round-off errors can accumulate 
in computations like the evaluation of high-degree polynomials. Such polynomials may have coeffcients with a large 
number of significant digits. Moreover, huge errors can accumulate if a polynomial is evaluated at a point whose 
value is a number with a large number of significant digits. 
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We can also solve inequality 7 using a combination of the Newton-Raphson (NR) method and Algorithm I. We 
call this method our Algorithm II. The main motivation here is to study the good behaviour of the function 
p(r) = r — q(r) and solve the inequality quickly even if we use a large starting value for r >> g(0). Our Algorithm 
II is as follows. Essentially, we use the function p(r) = r — q(r) and solve for p(r) — using NR iterations. We may 
start with a reasonably large value of r so that r > q(r) (or p(r) > 0). We may also start with a smaller value of r 
where p(r) is negative. We perform NR iterations until the iterate r does not change more than a small value say, 
10 -30 . The value of p(r) for the terminating value of r is observed to be negative in all our runs. Now we switch 
over to Algorithm I with this value of r as starting value. Note that Algorithm I always terminates when executed 
with such a starting value of r where p(r) is negative. Finally at termination, a sharp small bound r satisfying 
inequality 7 is obtained. The six tables show the results for six different polynomials. 

The function p(r) responds well to the NR method. As mentioned earlier, q'(r) is small (see Section EJ , 
making p'(r) close to 1 (note that p'(r) = 1 — q'{r)). Since p'{r) is close to 1, NR has quadratic convergence. We 
need only few NR iterations for each zero of each polynomial; we find that the value of r resulting after these small 
number of NR iterations is such that p(r) is negative or equivalently, r < q(r). As mentioned above, such a value 
of r is a suitable starting point for Algorithm I. 

Now we concentrate on the NR iterations in Algorithm II required to approximately compute the fixpoint of q. 
Let p = q{p) be the fixpoint. Suppose we start with a large starting value R for the NR iterate r. Note that in a 
few steps we would quickly get a value of the iterate r such that the first few significant bits of r match with the 
first few significant bits of p. This happens because the first step corrects the initial value R to R — p(R)/p'(R), 
which is nearly q(R) and very small compared to R. This happens since p{r) = r — q(r), p'(r) = 1 — q'(r) and 
q'(R) is small eventhough R may be quite large. So, in a few NR iterations r gets close enough to p to match p in 
a few significant bits (binary digits) . 

Suppose r and p have c > 1 identical (most) significant bits. Due to quadratic convergence property of NR 
iterations, the number of matching significant bits will increase in geometric progression with a ratio greater than 
unity. So, after k iterations, w k c bits would match, where w > 1. If we are interested in getting b bits of r match b 
bits of p, we would then need k steps where b = w k c. We assume without loss of generality that p = dx 2 E , where 
0.5 < d < 1 and E is any integer. Then, in k — 0{\og w b) iterations, we would get the iterate r close to p with 
relative error \r — p\/p < (2~ b x 2 E )/2 E ~ 1 = 2~ i - b ~ 1 K We summarize this in the following theorem 

Theorem 3.3 The fixpoint p — q(p) can be approximated within relative error 2~ b in O(logb) Newton-Raphson 
iterations, where b > 1 is an integer. 

3.4 Coping with errors in the input process 

Another aspect that needs mention is the possibility of small errors creeping in in the process of feeding inputs 
into our error computing algorithm. There are only two sets of inputs; the coefficients of the input polynomial 
and the approximate zeros. We can avoid input errors in the coefficients if each coefficient is integral. If we have 
decimal rational numbers as coefficients, the input of such numbers may suffer an error in the process of standard 
input in any programming language. So, we multiply all the coefficients with a suitable large power of 10 to make 
all coefficients integral. This does not change the polynomial and therefore the zeros remain unchanged by such 
multiplication of the coefficients of the input polynomial. Indeed, for this very reason, all examples (except example 
5) illustrated in this paper are about polynomials with integral coefficients. The approximate zeros in the input 
may however originate from any source or any numerical method. In the context of this paper, we require that 
approximate zeros be propagated accurately into our error finding Algorithms I and II; if there is any error in the 
input of approximate zeros, these errors must be accounted for in the final error bound computed for each zero. For 
the sake of simplicity, we may assume that such errors do not occur in the input of the approximate zeros; after all, 
once the intended approximate zeros are entered into variables of our error computing program, the values of the 
approximate zeros and all coefficients of the input polynomial are processed considering all errors that may occur 
during the computation of the final error bounds on each approximate zero. This is ensured by the floating-point 
filter and the real data type in LED A E3 ■ It is however possible and necessary to determine safe upper bounds 
in input processes for inputs from standard input instructions in different programming languages and systems and 
incorporate those error bounds suitably into the final error bound. A detailed study of the issue of errors in input 
processes for computations of safe boolean operations between polygons is done in [H]- Naturally, computations 
of errors resulting due to round-off errors in finite precision floating-point computations must also include due 
consideration of possible errors in input variables themselves; safe upper bounds on errors in each input variable 
must be provided and accounted for in the propagation of errors through the entire computation until outputs are 
generated. Such error analysis is reported in [15117], 

Example 3 

g(z) = 1000s 6 - 13016z 5 4- 59214z 4 - 107974z 3 + 61769z 2 - 997s + 4 
| Zeros at precision 7 | 



6 



tJounds using Algorithm 1 tor c = u.uuui 



Range zeros 
at precision 7 


Value Of 

9(0) 


Our bounds on Range 
zeros at precision 7 


Number of 
iterations 


1.993716 


2 .089829937337557E — 05 


2 . 0899 7i i J 2 '. )■> J ! !)-"• 


7 


4. 0080655 


3.64330, '2 l'>-,l!22:>E — 05 


;!. G 138 173 1573261E — 05 


14 


2. 997721 


1 . G fi 9 ( > ■' j ■ j S ! i ( i 1 - ■; Li it. , 


1 . 65977007537648E — 05 


7 


1 . 


7.32990974711188.E — 07 


7.329983!) ii):;. - ! /;." u7 


1 


.009422683 


1 .26751047717699E — 08 


1 . 267523 1 5228 1 76 S — 08 


1 


. 007075 149 


9.44998978980977B — 09 


9 . 45008428970767.E — 09 





Bounds using Algorithm II 



Range zeros 


Starting 
Value 


NR bounds on Range 


Number of 
NR + Algorithm I 


4.993716 


IE - 04 


2.09016807585603 - 05 


4+1 


4.0080655 


IE - 04 


3. 64417392553-1 -1 -IE - 05 


4+1 


2.997721 


0.01 


1.6599262404165E - 05 


4+1 


1.0 


0.01 


7.33067508160057E - 07 


4+1 


0.009422683 


0.01 


1.26764598079487E - 08 


4 + 1 


0.007075149 


IE - 10 


9.45098509041444E - 09 


3 + 1 



| Zeros at precision 16 ] 



Bounds using Algorithm I for c = 0.00000001 



at precision 16 


Value Of 
1(0) 


Our bounds on Range 
zeros at precision 16 


Number of 


4.99371584412958 


1.7990d3(iliil,s l()GE - 14 


1.7990436341745E - 14 


1 


4.00806551632572 


3.09602233507163E - 14 


3.09602236603186fi - 14 


1 


2.99772080758747 


1.38592484031255E - 14 


1.3859248541718E - 14 


1 


1.0 


5.77846972093955E - 16 


5.77846977872427E - 16 


1 


0.00942268285074248 


1.13201138948718E - 17 


1.1320114008073.E - 17 


1 


0.00707514910648869 


8.61633999170187E - 18 


8.61634007786531B - 18 


1 



Bounds using Algorithm II 



Range zeros 
at precision 16 


Starting 
value 


NR bounds on Range 
zeros at precision 16 


Number of 
NR + Algorithm I 


4.993715844129577 


IE - 13 


1.7990436341 746 IE - 14 


2 + 1 


4.008065516325719 


IE - 22 


3.09602236603224B - 14 


2 + 1 


2.997720807587473 


IE - 09 


1.38592485417189E - 14 


2 + 1 


1.0 


IE - 04 


5.77846977872429E - 16 


3 + 1 


0.009i22(.;S2Sr ) i}7121S.", 


IE - 04 


1.13201140080731E - 17 


3 + 1 


0.007075149106488685 


IE - 22 


8.61634007786535B - 18 


2 + 1 



Example 4 

g(z) = lOOOOOOz 4 + 223069s 3 - 41948404.Z 2 + 68883845s + 125362605 
| Zeros at precision 7 | 



Bounds using Algorithm I for e = 0.01 



Range zeros at 
precision 7 


Value at 

9(0) 


Our bounds on Range 
zeros at precision 7 


Number of 


4.00102 


1.14317714827099S - 03 


2.9420093174445E - 3 


95 


3.998911 


1.1428874099077E - 03 


2.93261103460139B - 3 


95 


-1.1 


2.11036257663918E - 07 


2.1314(>(i202 L\ir,r,HE - 7 


1 


-7.1229995 


5.45032569566617E - 07 


5.50482S!.>r,2(.i22S.l.E - 7 


1 



Bounds using Algorithm II 



Range zeros at 
precision 7 


Starting 
Value 


NR bounds on Range 
zeros at precision 7 


Number of 
NR + Algorithm I 


4.00102 


IE - 9 


2.93304493819802fi - 3 


9 + 1 


3.998911 


IE - 1 


2.9328822051G119B - 3 


9 + 1 


-1.1 


IE — 5 


2.11057117059151E' - 7 


3 + 1 


-7.1229995 


IE — 5 


5.45087277340316E - 7 


3 + 1 



7 



Bound, using Algorithm I for c = 0.00000001 



Ti.aii'-^e zeros at 


Value at 

9(0) 


Our bounds on Range 
zeros at precision 16 


Number of 
iterations 


4.001019657681204 


1.47961652695472E - 12 


1.47961G.V117r,0HS.E - 12 


1 


3.998910834038721 


1.47927388424087E - 12 


1.47927389903361E - 12 


1 


-1.100000011624190 


2.3844189558305E - 16 


2.34844191906748.E - 16 


1 


-7.122999480095735 


6.67094208084199E - 16 


6.67094214755143 - 16 


1 



Bounds using Algorithm II 



Range zeros at 
precision 16 


Value 


NR bounds on Range 


Number of 
NR + Algorithm I 


4.001019657681204 


IE - 11 


1.47961654279053B - 12 


3 + 1 


3.998910834038721 


IE - 22 


1.47927390007277B - 12 


2 + 1 


-1.100000011624190 


IE - 11 


2.34844191906748.E - 16 


2 + 1 


-7.122999480095735 


IE — 5 


6.67094214755114.E - 16 


3 + 1 



Example 5: 

g(z) = 8.7029s 9 - 167s S + 463.33s 7 + 1126. Is 6 + 76.241s 5 - 7.0508s 4 - 4085.4s 3 - 1036.1s 2 - 99.729s - 54.649 



| Zeros at precision 7 | 
Bounds using Algorithm I for e = 0.0001 



Actual 


Range zeros 
at 

precision 7 


Value of 

9(0) 


Our bounds on 

Range zeros 
at precision 7 


Number 
of 

iterations 


15.0911133 


15.09111 


9.076342589 
63854E - 06 


9.07725022 
389751E - 06 


1 


5.60181709 


5.601817 


5.206365229 
73372E - 06 


5.20688586 
625669B - 06 


1 


1.471636 


1.471636 


7.629781518 


7.6305444 


1 






41216E - 07 


96564E - 07 




0.029598 + 
iO. 204796 


0.02959805 
+ i0. 2047964 


8.494188735 
02682E - 08 


8.49503815 
390032E - 08 


1 


0.029598- 
iO. 204796 


0.02959805 
-iO. 2047964 


8.494188735 
02682E - 08 


8.49503815 
390032E - 08 


1 


-0.3115666 


-0.3115659 


1.57728567 
206155E - 07 


1.5774434 
0062876S - 07 


1 


-0.396041 
3425333 


-0.396041 
+ 342533 


7.26872135 
734945E - 07 


7.26944822 
958519E - 07 


1 


-0.396041 
-il. 3425333 


-0.3960412 
-il. 342533 


7.26872135 
734945E - 07 


7.26944822 
958519E - 07 


1 


-1.9311091 


-1.931109 


1.03460887 
418129E - 06 


1.034712335 
06871E - 06 


1 



Bounds using Algorithm II 



Range zeros 
at 

precision 7 


Starting 
Value 


Our bounds on 
Range zeros^ 


Number 
of 

iVff+Algorithm I 


15.09111 


IE - 06 


9.07733851 
50S154E - 06 


3 + 1 


5.601817 


IE ~ 06 


5.2069566 
1442314S - 06 


3 + 1 


1.471636 


IE - 11 


7.63058385 
852901E - 07 


3 + 1 


0.02959805 
+ i0. 2047964 


IE - 11 


8.49504782 
584932E - 08 


3 + 1 


0.02959805 
-iO. 2047964 


IE - 02 


8.49504782 
584932E - 08 


4+ 1 


-0.3115659 


IE - 05 


1.57744695 
490493E - 07 


3 + 1 


-0.396041 
342533 


IE - 02 


7.2694861 
2010097S - 07 


4+ 1 


-0.3960412 
-il. 342533 


IE - 03 


7.2694861 
2010097S - 07 


4+ 1 


-1.931109 


IE ~ 03 


1.03471867 
334759E - 06 


4+ 1 



8 



I Z^eros ai precision 1 b | 
Bounds using Algorithm I for e = 0.00000001 



at precision 16 


Value of 

9(0) 


Our bounds on 


Number 
of 


15.09111310576714 


5.53750340 
672293E — 15 


5.53750346 
209812E — 15 


1 


5.601817329197456 


3.14720692 
264054E — 15 


3.14720695 
411265E — 15 


1 


1.471636442949784 


5.8533007 
612049E — 16 


5.85330076 
12049 E — 16 


1 


0.02959805472609926 
+ 7 . 2 1 7 9 6 1 1 5 9 305099 


6.97248888 
506532E — 17 


6.972488954 
79027E — 17 


1 


0.02959805472609926 
— 70 . 20 1796 1 159305099 


506532E — 17 


79027E — 17 


1 


-0.3115658643893903 


1.461254491 
76875E — 16 


1.46125450 
63813E — 16 


1 


-0.3960411556175865 
+ 7 1.34253271 7206571 


5.65575203 
517524E - 16 


5.655752091 
73281E - 16 


1 


-0.3960411556175865 
-71. 342532717206571 


5.65575203 
517524E - 16 


5.655752091 
73281E - 16 


1 


-1.931108848212619 


7.35055936 
523337E - 16 


7.350559365 
23337E - 16 


1 



Bounds using Algorithm II 



Range zeros 
at precision 16 


Starting 
Value 


Our bounds 
at precision 16 


Number 
of 

NR + Algorithm I 


15.091113105767141 


IE - 14 


5.53750346209816E - 15 


2 + 1 


5.6018173291974556 


IE — 7 


3.14720695411267E - 15 


3 + 1 


1.471636442949784 


0.0001 


5.85330076120492E - 16 


3 + 1 


0.02959805472609926 
+ ■70. 2047964459305099 


IB - 22 


6.97248895479027E - 17 


2 + 1 


0.02959805472609926 
-70. 2047964459305099 


IE - 14 


6.97248895479027E - 17 


2 4- 1 


-0.3115658643893903 


IE — 7 


1.46125450638131E - 16 


3 + 1 


-0.3960411556175865 
+ 71. 342532717206571 


0.0001 




3 + 1 


-0.3960411556175865 
-7 1.34253271 7206571 


IB - 14 


5.65575209173284E - 16 


2 + 1 


-1.931108848212619 


IE - 22 


7.3505593652334E - 16 


2 + 1 



Example 6: 

g(z) = z 20 + 10 12 z 14 + ; B + 1 
r Zero, at precision 7 | 
Bounds using Algorithm I for e = 0.0001 



Range zeros 
at precision 7 


Value 
of 

9(0) 


Our bounds on 
Range zeros 


Number 
of 

iterations 




■i.(iO:i02117122171E - 7 


4.66410117625052E — 7 


1 


86.60254 - 750 


4.66363 18127691 IE - 7 


4.66410117625041E — 7 


1 


0.1354659 + iO. 03091968 


6.1585340861066E - 8 


6.15914993951521E - 8 


1 


0.1354659 - iO. 03091968 


6.1585340861066E - 8 


6.15914993951521E - 8 


1 


0.1086348 + iO. 0866332 


6.15859858647888E - 8 


6.15921444633753E - 8 


1 


0.1086348 - zO.0866332 


6.15859858647888E - 8 


6.15921 1 1 1633753E - 8 


1 


0.06028841 + 7.0. 1251894 


6.15851872()()21,siE - 8 


6.15913457189681E - 8 


1 


0.06028841 - 70.1251894 


6.15851872002181E - 8 


6.15913457189681E - 8 


1 


1.666667E - 29 + ilOO 


4.6636348316034E - 7 


4.66410119508656E - 7 


1 


1.666667E - 29 - 7100 


4.6636348316034E - 7 


4.66410119508656E — 7 


1 


-5.140612E - 7 + 10.1389 195 


6.15855211502818E - 8 


6.15916797023968E - 8 


1 


-5.140612E - 7 - 70.13891!).-, 


6.15855211502818E - 8 


6.15916797023968E - 8 


1 


-0.0602879 . (0.125189 


6.15858807166052E - 8 


6.15920393046768E - 8 


1 


-0.0602879 - 70.125189 


6.15858807166052E - 8 


6.15920393046768E - 8 


1 


-0.1086354 + 10.0866340 


6.15850868338623E - 8 


6.15912 153 125 157E - 8 


1 


-0.1086354 - 70.0866340 


6.15850868338623E - 8 


6.15912 153 125 156E - 8 


1 


-0.1354657 + 70.03091868 


6.15857290119851E - 8 


6.15918875848863E - 8 


1 


-0.1354657 - 70.03091868 


6.15857290119851E - 8 


6.159l8,s75,Mfsfs03E - 8 


1 


-86.60254 + 750 


4.66363481988736E - 7 


4.66410118336934E - 7 


1 


-86.60254 - 750 


4.66363481988736E - 7 


4.66410118336934E - 7 


1 



9 



riounds using Algoritnm 11 



Range zeros 


Starting 
Value 


NR bounds on 
Range zeros 


Number 
of 

iterations 
iV.R+ Algorithm I 


86-60254 + i50 


IE - 10 


4.66410190829221.E - 7 


2 + 1 


86.60254 - i50 


IE - 06 


4.66410190829219.E - 7 


2 + 1 


0.1354659 + i0.03091968 


IE - 05 


6.15920100912469.E - 8 


3 + 1 


0.1354659 - iO. 03091968 


IE - 06 


6.15920095777814.E - 8 


2 + 1 


0.1086348 + iO. 086633225 


IE - 10 


6.15926551678678E - 8 


2 + 1 


0.1086348 - iO. 086633225 


IE - 05 


6.15926551700718.E - 8 


3 + 1 


0.06028841 + zO.1251894 


IE - 05 


6.1591856-1125151.E - 8 


3 + 1 


0.06028841 - iO. 1251894 


IE - 03 


6.15918564125148.E - 8 


3 + 1 


1.666667S - 29 + ilOO 


IE - 04 


4.66410192449928.E - 7 


2 + 1 


1.666667E - 29 - ilOO 


IE - 03 


4.66410192712841.E - 7 


3 + 1 


-5.140612E - 7 -)- 10.1389-1955 


IE - 04 


6.15921904014471.E - 8 


3 + 1 


-5.140612E - 7 - 10.1389-1955 


IE - 06 


6.159218988798E - 8 


2 + 1 


-0.0602.x 79 ■ (.0.125189 


IE — 03 


6. 159255mm! Hi 11 (>£ - 8 


3 + 1 


-0.0602879 - iO. 125189 


IS - 06 


6.1592519i9C1725£ - 8 


2 + 1 


-0.1086354 + iO. 08663403 


IE - 04 


6.15917560344614E - 8 


3 + 1 


-0.1086354 + i0.08663403 


IE - 10 


6.15917560322547E - 8 


2 + 1 


-0.1354657 + 10.03091868 


IE — 06 


6.1592397773!.H)59.E - 8 


2 + 1 


-0.1354657 - iO. 03091868 


IE - 10 


[6.15923982851674E - 8 


2 + 1 


-86.60254 + i50 


IE - 04 


4.66410191278206.E - 7 


2 + 1 


-86.60254 - i50 


IE - 10 


4.66410191541114E - 7 


2 + 1 



| Zeros at precision 16 | 
Bounds using Algorithm I for e = 0.00000001 



Range zeros 
at precision 16 


Value of 

9(0) 


Our bounds 
at precision 16 


Number 
of 


86.6025-10378-1 1386 + i50 


2.77298663308076E - 15 


2.77298666081066E - 15 


1 


86.60254037844386 - i50 


2.77298663308073E - 15 


2.77298666081066E - 15 


1 


0. 135465908523133 
+i0. 03091968468442907 


8.26333010777928E - 17 


8.26333019041465E - 17 


1 


0.135465908523133 
-iO. 03091968468442907 


8.26333010778121E - 17 


8.26333019041465E - 17 


1 


0.1086348117116894 
+i0. 0866332513523468 


8.26340929550181E - 17 


8.26340937813606E - 17 


1 


0.1086348117116894 
-iO. 08663325 13523468 


8.26340929550181E - 17 


8.26340937813606E - 17 


1 


+i0. 125189441295516 




8.26331063087529E - 17 


1 


0.06028841320891334 
-iO. 125189441295516 


8.26331054824203E - 17 


8.26331063087529E - 17 


1 


1.666666666666667E - 29 + £100 


2.77298663308073E - 15 


2.77298666081066E - 15 


1 


1.66666666666667E - 29 - ilOO 


2.77298663308073E - 15 


2.77298666081066E - 15 


1 


-5.140611950289672E - 7 
+i0. 1389495494401665 


8.26335449713578E - 17 


8.26335457976949E - 17 


1 


-5.14061950289672E - 7 
-iO. 1389495494401665 


8.26335449713578E - 17 


8.26335457976949E - 17 


1 


-0.06028748690264759 
+i0. 1251889952099292 


8.26339844245486E - 171 


8.26339852508899E - 17 


1 


-0.06028748690264759 
-iO. 1251889952099292 


8.26339844245486E - 17 


8.26339852508899E - 17 


1 


-0.1086354527355147 
+i0. 08663402895368696 


8.26329969321128E - 17 


8.26329977584442E - 17 


1 


-0.1086354527355147 
- iO . 08663402895368696 


8.26329969321128E - 17 


8.26329977584442E - 17 


1 


-0.135465797443785 
+10.03091868233921389 


8.26337888538945E - 17 


8.2633789680234E - 17 


1 


-0.1354657443785 
-iO. 0309 186823392 1389 


8.26337888538945E - 17 


[8.2633789680234E - 17 


1 


-86.60254037844386 + i50 


2.77298663308073E - 15 


2.77298666081066E - 15 


1 


-86.60254037844386 + i50 


2.77298663308073E - 15 


2.77298666081066E - 15 


1 



10 



.Bounds using Algorithm 11 



at precision 16 


Starting 
Value 


NR bounds 
at precision 16 


Number 
of 

iterations 
JVR + Algorithm I 


86.60254037844386 + i50 


IE - 15 


2.77298666081066E - 15 


2 + 1 


86.6025 iu:;7S i i:-J86 - i50 


IB - 10 


2.77298666081066-E - 15 


2 + 1 


0.135465908523133 
+ i0. 03091968468442907 


IE - 07 


8.26333018961745.E - 17 


2 + 1 


0.135465908523133 
-iO. 0300101,8 108 1-12907 


IE - 15 


8.26333019041475B - 17 


1 + 1 


0.1086348117116894 
+i0. 0866332513523468 


IE — 10 


8.26340937813615B - 17 


2 + 1 


0.1086348117116894 
-iO. 0866332513523468 


IE - 16 


8.26340937813615.E - 17 


1 + 1 


0.06028841320801334 
+i0. 125189441295516 


IE — 10 


8.26331063087538E - 17 


2 + 1 


0.06028841320891334 
~i0. 125180441295516 


IE — 07 


8.26331063007808B - 17 


2 + 1 


1.666666666666667.E - 29 + ilOO 


IE - 15 


2.77298666081066S - 15 


2 + 1 


1.666660i>i>i>i>i>i>i>i>7E - 29 - ilOO 


IE — 10 


2.77298666081066-E - 15 


2 + 1 


-5.140611950289672.E - 7 
+i0. 1389495494401665 


IE - 10 


8.26335457976958.E - 17 


2 + 1 


-5.140(il0.",i}280(i72£; - 7 


IE — 07 


8.26335 157fs0722GE - 17 


2 + 1 


-i0. 1389495494401665 








-0.06028748690264759 
+ i0. 1251889952099292 


IE - 16 


8.26339852508908.E - 17 


1 + 1 


-0.06028748690264759 
-iO. 1251889952090292 


IE - 07 


8.26339852429179.E - 17 


2 + 1 


0. 1 (i>-.'. ir,27:i55147 


IE - 15 


8.26329977584451E - 17 


1 + 1 


+ i0. 08663402895368696 








-0.1086354527355147 
-iO. 08663402895368696 


IE - 16 


8.26329977584451.E - 17 


1 + 1 


-0.135465797443785 


IE - 10 


8.26337896802349.E - 17 


2 + 1 


+i0. 03091868233921389 








-0.135465797443785 


IE - 15 


8.26337 j.!! a: ; 7 


1 + 1 


-iO. 03001868233921389 








-86-602T,. UK!7S 1 K5S6 + i50 


IE - 16 


2.77208(,(,(,l)81(K,6A7 - 15 


2 + 1 


-86.6025 10:i7fs 11:180 - i50 


IE - 07 


2.7720> ::-.;::,.,;/•; - 15 


2 + 1 



4 Results and observations 

We tabulate the bounds computed for examples 1 through 6 using our Algorithms I and II based on Theorem 12 .11 
We also tabulate the number of iterations required by Algorithms I and II. For Algorithm II, the Newton-Raphson 
iterations are followed by steps of Algorithm I. So, the numbers of iterations of both these stages are depicted in 
the tables. Note that the value of q(0) for each approximate zero is very close to (and slightly lesser than) the error 
bound, as argued in Section l3~Tl 

For examples 1 and 2 the actual zeros are known. Both these polynomials (see have distinct and well 
separated zeros. Approximate zeros as computed by the ZERPOL algorithm ^H] are used in p<3; Smith computes 
error bounds for these zeros using his method (as in JHj) based on Gerschgorin's theorems. Error bounds computed 
by our Algorithm I using Theorem 12.11 for these ZERPOL zeros are comparatively inferior to those of [Ej- F° r 
example 1, our bounds are comparable to those of Smith; the bounds on errors in the last two zeros are very 
close. For example 2, our bounds for the seventh and eighth zeros are in fact better than those obtained in |19j . 
The bounds for the other eight zeros are almost of the same order in our case whereas Smith's bound for these 
eight zeros are much better. However, our method is scalable. We show that our bounds are sharper for more 
accurate approximations. As we increase the number of significant digits in the approximate zeros computed using 
the Range software [T|, we find that our error bounds too improve as expected. We use the Range software for 
computing approximate zeros for all examples at precisions of 7 and 16 decimal significant digits, respectively. 

Our method works well for polynomials with close zeros as shown in examples 3 and 4. The close zeros in 
example 3 are smaller than the rest of the zeros. Error bounds computed using our Algorithms I and II for zeros 
computed using Range at the two precision values of 7 and 16 are presented. The bounds in example 3 for small 
and close zeros are much sharper than the bounds for the other zeros because the numerator l(r) of q(r) becomes 
smaller for smaller zeros. In contrast, we see a coarser bound for close zeros with larger magnitudes in example 4, 
even though they are as close as the close zeros in example 3. 

The bounds for the two close zeros (first two zeros) in example 4 at precision 7 are larger than the separation 
between the zeros. The two circles with radii equal to these two computed bounds and centred at the respective 
approximate zeros enclose both the approximate zeros. This is guaranteed by Rouche's theorem because the 
condition r > q(r) holds for the bound r for these two zeros. In other words, the two approximate zeros and the 
two corresponding actual zeros lie in both the above mentioned circles. Our bounds for the more separated zeros 
in example 4 are better that those for the closer zeros. Increasing the accuracy of approximation helps in getting 
sharp bounds isolating each zero. This is observed even for examples 3 and 4 at precision 16. 

Example 5 is from and example 6 is from ^2]- These are relatively higher degree polynomials and our 
Algorithms I and II show consistent and good bounds. We make an crucial observation over all examples; we note 
that the relative error in approximate zeros as computed by our Algorithms I and II are roughly of the same order 
of magnitude over all zeros for each example. 
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Posteriori error bounds on approximate zeros of univariate polynomials may be used in geometric computations 
with high degree curves and surfaces as well as in various scientific, engineering and computational mathematics 
applications. Our method for computing error bounds can be used in algorithms for computing approximate zeros 
for polynomials to desired accuracies; for instance, a Newton-Raphson based algorithm for computing approximate 
zeros may repeatedly compute error bounds in each iteration using our method to decide whether approximate zeros 
of desired accuracy have already been computed. Use of high precision in zeros' computations can help generating 
as close approximations to actual zeros as one wishes. Krishnan et al. [2] use Durand-Kerner iterations for fast 
convergence to all zeros of a polynomial simultaneously; at each step they check the quality of approximation 
achieved up to that step using Smith's method JSj for computing upper bounds on errors in the approximations 
of zeros obtained. Our method of computing similar error bounds as shown in this paper can be used in the place 
of Smith's method. 

We have used requisite high precision in our algorithms for decision making steps involving inequalities and in 
the evaluation of polynomials of high degree. Our method works for polynomials with close clustered zeros; the 
problem in such cases is that we might require larger numbers of significant digits in the approximate zeros for 
computing good bounds. With coarser approximations, we might get poorer bounds (see example 4). However, 
for the case of close roots with small magnitudes, as in example 3, we may get sharp bounds even with coarse 
approximations of zeros. We can also extend our method for handling the case of polynomials whose zeros have 
multiplicity greater than unity. The bounds obtained for such cases are likely to be coarse and we feel that only 
very precise approximations of zeros would give sharper bounds. 
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